Topological surface states scattering in Antimony 
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In this work we study topologically protected states on the Sb(lll) surface by using ab-initio 
transport theory. In presence of a strong surface perturbation we obtain standing wave states 
resulting from the superposition of spin polarized surface states. By Fourier analysis, we identify 
the underlying scattering processes and the spin texture. We find evidence of resonant transmission 
across surface barriers at quantum well states energies and evaluate their life-times. Our results are 
in excellent agreement with experimental findings. We also show that despite the presence of a step 
edge along a different high symmetry direction, not yet probed experimentally, the surface states 
exhibit unperturbed transmission around the Fermi energy for states with near normal incidence. 

PACS numbers: 



In the last few years topological insulators (TIs) have 
become an intensely studied field of condensed matter 
physics 0, 0] . The unique metallic states at the surface of 
a TI can be used both as a table-top playground to prove 
fundamental concepts developed for particle physics and 
also as new materials platform for intriguing quantum 
applications in both spin-electronics and quantum com- 
puting ■ The first three-dimensional TI to be exper- 
imentally realized was a bismuth-antimony alloy [6| , fol- 
lowing a theoretical prediction by Fu and Kane [7j. The 
topological nature of the alloy is inherited from the par- 
ent element antimony, which has a non-trivial principal 
topological invariant i/q, whereas bismuth has a trivial 
invariant Q. 

Although Sb itself is a semimetal, its (111) surface 
hosts two spin-polarized bands, which extend around the 
Fermi energy, E-p. These form a single distorted Dirac 
cone where the lower bands are lifted upwards. The elec- 
tronic properties of Sb(lll) thin films have been investi- 
gated both theoretically and experimentally In a 
beautiful experiment, Seo et al. demonstrated that the 
topological surface states are extraordinarily insensitive 
to the presence of surface barriers fl3l ] . They probed 
the extended nature of Sb(lll) surface states by using 
a scanning tunneling microscope (STM) and found that 
these transmit across surface atomic steps with a high 
probability. Their analysis of the standing-wave states 
on surface terraces revealed the novel chiral spin texture 
of the two surface states, consistent with earlier Angle- 
resolvedphotoemission spectroscopy (ARPES) measure- 
ments 
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In this Letter we seek to theoretically recreate the 
above mentioned experiment by using ab-initio transport 
theory, and show that we can reproduce the formation of 
quantum well states and their life-times, as well as the 
wave lengths and phase shifts of the scattering states. 
Thereby we demonstrate that by first principles calcula- 
tions one can describe the correct scattering properties of 
such topologically protected surface states. In addition 
of comparing our results favorably to the experiments, 



we predict the scattering properties of these states in 
presence of a surface perturbation along a direction or- 
thogonal to the one probed experimentally. 

Our first-principles electronic structure calculations 
are performed with density functional theory (DFT) 
using the local density approximation (LDA) to the 
exchange-correlation functional. We employed the 
SIESTA package, which implements a linear combination 
of atomic orbitals basis set flBj ]. Spin-orbit interaction, 
essential to describe the surface states, has been included 
via the on-site approximation flil ]. The transport prop- 
erties are then calculated by using SMEAGOL, which 
combines the non-equilibrium Green's function (NEGF) 
method with DFT (T^-El- Our order- N implementation 
of the code allows us to treat large systems [2(|. We use 
a double-C polarized (DZP) basis set, with a cutoff en- 
ergy of 300 Ry for the real space mesh. We have carefully 
checked convergence of our results with respect to all the 
parameters used. 

Antimony crystallizes in a rhombohedral structure 
[space group D'^ d (R3m)] with two atoms per unit cell. 
An alternate way to represent its structure is in an hexag- 
onal setting with the unit cell comprising six atoms. This 
representation is particularly useful to construct two- 
dimensional slabs, which arc made of Sb bilayers, as 
shown in Fig. [TJa). The inter-bilayer coupling is weak 
and it is possible to create surface steps which are a sin- 
gle bilayer high [l3l [Hf . The bulk structure was relaxed 
using the Vienna ab-initio simulation package (VASP) 



2l| , until the forces were less than 0.01 eV/A. 



We begin by calculating the surface band structure of 
six and twelve bilayers thick slabs of Sb [Figs. [Ub) and 
(c)], by using a 10 x 10 in-plane fc-point grid. The dis- 
torted Dirac cone at T is gapless indicating minimal inter- 
action between the top and bottom surfaces of the slab. 
The surface band structure matches well previous ab- 
initio calculations 

0,113, El- We find the Dirac point at 
an energy of about 210 meV below Ef for six and twelve 
bilayer slabs. In order to simulate the ARPES spectrum 
of an infinitely thick slab we perform a SMEAGOL calcu- 
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FIG. 1: (Color online) (a) Antimony unit cell in the hexag- 
onal setting. The atoms form a bilayer structure with the 
intralayer distance as 1.51 A and the interlayer distance as 
2.25 A. Band structure for (b) six bilayers and (c) twelve bi- 
layers thick slabs along K — T — M directions, (d) Simulated 
ARPES from a semi-infinite slab. The distorted Dirac cone is 
found to be comprised of spin-polarized surface bands. Here 
warmer colors represent higher PDOS (red represents largest 
values, blue lowest ones, with the color scale in between be- 
ing linear). Spin-resolved ARPES along (e) f - M and (f) 
F — K directions showing the opposite spins of the two surface 
bands along the directions indicated by arrows in the inset of 
the figures. In this case, red and blue colors indicate up and 
down spins, respectively. 



lation for the 6 bilayer slab, where we attach semi-infinite 
Sb electrodes at the bottom layer via self-energies. The 
ARPES spectrum is then obtained by calculating the pro- 
jected density of states on the surface atoms, and the 
result is shown in Fig. (Ud). The spin- resolved ARPES 
[Fig.QJc) and (f)] shows that the two surface bands carry 
opposite spin and exhibit the characteristic spin texture 
associated with topologically non-trivial materials. Fur- 
thermore, it can be seen that the surface states are more 
pronounced close to the T point. Once the two bands 
turn downwards from their maximum point, they are less 
localized on the surface due to their hybridization with 
the bulk bands. This matches with ARPES experiments, 
where at M, no surface states are found 



1J. Thus, in 



both scattering and ARPES experiments, one would ex- 
pect the dominant contributions to come from an area 
around the center of the BZ, with a radius of about one 
third of the length of BZ along T — M direction. 

Next we simulate a step perturbation on the surface 
for two different directions: in the first case the step runs 
along the y-direction [see Fig. HJb)], so that the trans- 
port is along the ^-direction, which is the same as in the 
experiment in Rcf. [131 ] : and as a second orientation we 
choose the orthogonal direction (step running along x), 



to evaluate the effect of the orientation on the scattering. 
We have relaxed the step geometry for a smaller scatter- 
ing region, but the atomic displacements were only mi- 
nor, and henceforth we use the unrelaxed step configura- 
tion. A single bilayer high step is created over a length of 
120 A. The adjacent flat region extends over 270 A. The 
setup consists of a 13 bilayer-thick region with a short 12 
bilayer-thick region on the left and a longer one on the 
right, as shown schematically in Fig. [2f a) ( for the actual 
setup see Supplemental Material (SM) [23J). We attach 
semi-infinite leads on left-hand and right-hand sides of 
the scattering region, by means of the self-energies calcu- 
lated by SMEAGOL. This enables us to simulate isolated 
scatterers, thus avoiding the use of periodic boundary 
conditions along the transport direction [2fjj . 

The total projected density of states (PDOS) is ob- 
tained by integrating over all fc^-points perpendicular to 
the transport direction 



PDOS to tal(^) 



P~DOS k± (E) dkj 



(1) 



Analogously the total transmission is given by 
TtotaK-E 1 ) = J k± T k± (E)dk A _. We note that in the 
first orientation of the step k± runs along the T — K 
direction and the transport direction is parallel to T — M 
in reciprocal space, whereas in the second orientation k± 
runs along T — M and the transport direction is parallel 
to f — K . While we find that 3 /c^-points are sufficient 
for obtaining a converged self-consistent potential, we 
need many more fc^-points for accurately integrating 
PDOStotai for a given potential, where we therefore use 
200 fc^-points. 

The PDOSfc^ of the atoms on the top surface for 
k± = is shown in Fig. [2jc). The quantum well states 
formed by quantization of the energy levels in the step 
region are clearly visible, and extend over the energy win- 
dow in which the two topological surface bands exist. 
The PDOS on the adjacent flat region shows oscillatory 
behavior typical of one-dimensional scattering barriers 
[20| , and it has phase shifts at energies corresponding to 
the allowed energy states in the step region. At those en- 
ergies wc also find resonant transmission across quantum 
well states, visible as peaks in the transmission curve in 
FigSd). This matches the experimental observation of 
resonant tunneling through the surface barrier at those 
energies. It indicates a remarkably long phase coherence 
length of hundreds of angstroms, which is due to the ex- 
tended nature of topological surface states [T3|. Further, 
the change of phase when the states are reflected from the 
barrier is nearly zero. Over the entire energy range there 
is clearly a rather large amount of scattering caused by 
the step and the transmission drops significantly below 
the value in absence of step. Wc have verified that the 
single bilayer step at the top surface perturbs the bottom 
layer minimally and there is only small coupling between 



the two surfaces even in the presence of the surface step 
0. 

The integrated PDOS to tai is shown in Fig. [2jc). The 
main features corresponding to the quantized energy lev- 
els can still be identified at almost the same energies 
found for k±_ = 0, but are broadened and less pronounced. 
The broadening of energy levels for k± = is 8-12 meV, 
while for PDOS to tai it increases by a factor of about 
three. Hence, the lifetimes, which are inversely propor- 
tional to the broadening, are in the same range as those 
obtained experimentally (Tj|. This shows that the scat- 
tering properties of the step are well-reproduced in our 
calculations. We note that quantitatively the results for 
a 12 bilayer slab differ somewhat from the ones one would 
obtain for a semi-infinite surface due to the different sur- 
face band structures [Fig. [Tfc) vs. Hid)]. Firstly, the 
linear dispersion region of the surface states around T 
is more extended in energy for the semi-infinite surface, 
and this results in quantum well states spacing remain- 
ing constant, whereas for the 12 bilayer slab the spacing 
reduces with increasing energy. Moreover for 12 bilayers 
we find band edges at -185 meV and -108 meV, which 
are absent in the semi-infinite system, and which give 
strong contributions to the 12 bilayer PDOS. Due to the 
enhanced PDOS at these band edges the phase shifts cor- 
responding to resonant tunneling are not as clearly pro- 
nounced for PDOStotai as for k± = 0. Finally, whereas 
for the semi-infinite slab there are only two discrete sur- 
face bands around Ep, and all the other bands are dif- 
fuse, for the 12 bilayer slab clearly the number of distinct 
bands is larger, which as we will show leads to many more 
features in the scattering. These lead to discrete short- 
wavelength scattering processes which would be absent 
for a semi-infinite slab. 

In order to analyze the scattering mechanism and 
probe the chiral spin texture of the surface bands, we 
evaluate the Fourier transform (FT) of PDOS to tai along 
the flat region on the right side of the step, and the re- 
sult is shown in Fig. Hff). For a given energy and k± 
(k± is conserved during scattering) there are standing 
waves with a wave vector q for every pair of scattering 
states in the Sb electrodes with opposite group velocities. 
For each such pair of states with indices i and j we get 
the norm of scattering vector = | fcj — kj | , at which 
therefore the FT has enhanced amplitude. In Fig. (5Jf) 
we find a rather large amount of noise at low q, partly 
due to the fact that we integrate PDOStotai over a finite 
number of ^-points, nevertheless we can identify two 
prominent features that are preserved for increasing q. 
The first is an enhanced amplitude starting at about - 
210 meV and then increasing linearly with q with a slope 
of 1.1 eVA, and the second is an equivalent enhancement 
with the similar slope starting at -110 meV. These two 
prominent scattering wave-vectors are also found in the 
experiments of Seo et al., where they are labeled qs and 
qA, respectively, qs is attributed to scattering between 
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FIG. 2: (Color online) (a) A schematic of the setup showing 
one bilayer high surface perturbation on a twelve bilayer slab, 
(b) Top view of the slab; surface step runs along y for the first 
setup and along x for the second, (c) PDOS for surface atoms 
with a single step adjacent to a flat region for the first setup 
at k± = 0. The quantization of the energy levels in the step is 
clearly seen, with accompanying phase shifts in the adjacent 
flat region, (d) Transmission at k± = with and without the 
step indicating finite scattering due to the step. Average of 
transmission over all fcx-points is also shown, (e) PDOS for 
surface atoms averaged over all k±. (e) Fourier transform of 
PDOS data in (d) over the flat region reveals the different al- 
lowed scattering wave vectors. The most prominent features, 
qA and qb, are nearly linear with slopes equal to 1.1 eVA. 



the surface states close to T having the same spin, but 
opposite momentum direction. In this case the scatter- 
ing state momenta are not equal, unlike in conventional 
scattering, because of the unique spin texture of the sur- 
face states resulting in an asymmetric band structure for 
a given spin. The other scattering wave- vector (g^) is 
attributed to scattering between neighboring hole pock- 
ets away from T. The calculated slopes of qA and qs are 
remarkably close to the value of 1.2 eVA found experi- 
mentally [13j . 

An STM experiment measures the total scattering in 
all the reciprocal space directions. In order to establish 
which k± -points give rise to the qA and qs features in the 
average, we decompose the allowed scattering processes 
along different k± directions. To illustrate the general 
scattering mechanism, in Fig. (3Jc) we show the Fourier 
transform for an arbitrary &;^-point, and the bands of the 
Sb electrodes along the transport direction for the same 
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k±. As an example, we consider the energy at -150 meV, 
where in the band structure we find 2 bands crossing 
for positive fc, at k\ and fc 2 . In the Fourier transform 
we find 3 scattering vectors: q\ = k 2 — fci, q 2 = 2fci 
(obtained by k\ scattering to — fci), and similarly = 
2k 2 . The amplitude is highest towards the band edge, 
where the PDOS is maximal. Analyzing the FT for all 
fc^-points, we find that most such features found for a 
single k±_ disappear when averaging, except for qA and 
qs- There are two key factors, which decide which q 
vectors dominate. Firstly, for a given k± there needs to 
be a high scattering between the initial (A;,) and final (kj) 
states. This is the case if the spins in the two states are 
aligned and if the surface-PDOS for them is large. This, 
therefore, excludes bulk states and states with opposite 
spins. Secondly, for the features to remain prominent 
when integrated over k± , it is necessary that there is an 
extended region in the BZ where these features are found 
more or less unchanged. This is true when the change of 
band structure along k± is small, which is the case close 
to a maximum or a minimum, i.e. when dE/dk± = 0, 
and ideally when the band curvature is small (d 2 E/dk\ 
is small) . We have verified that the above two conditions 
are indeed satisfied for q^ and qs- From the 12 layer slab 
band structure along F — K [Fig. Hfc)] we identify two 
such features: the first is the minimum at about k± = F, 
and the second is the maximum at k±b = 0.167T. In Figs. 
[Ha) and (b) we show the scattering wave vectors for these 
two k± . We clearly identify the features leading to qA and 
qB- The remaining features in the FT disappear under 
averaging, since for these the aforementioned conditions 
are not satisfied. 

Finally, we create a step on the Sb(lll) surface along 
the rr-direction. In this case the slab consists of six bilay- 
ers, since simulating a 12 bilayer scattering region along 
this direction is beyond our computational resources (see 
SM for a description of the setup [3]). This direction 
is promising since there is an energy window [-60 meV 
to 20 meV in Fig. [IJb)] over which only a single spin- 
polarized surface state exists for k± = (note that in 
this case the transport direction is parallel to the f — K 
line in reciprocal space), which is reminiscent of the pro- 
totypical TI Bi2Sc3. The PDOSfc^ of the surface atoms 
for k± = is plotted in Fig. Ufa). In the energy range of 
a single spin-polarized state there is no scattering, which 
is the hallmark of a topological surface state. From - 
170 meV to 70 meV quantum well states are found as 
a result of superposition between the two surface bands, 
with a mechanism analogous to that in the first step ori- 
entation. Fig. |U[b) shows the transmission at k± = 
with and without the surface step. Remarkably, there 
is a perfect transmission in this energy window, despite 
the presence of the strong surface perturbation in the 
form of an extended single bilayer high surface step. This 
can be explained by invoking the general principle that 
disorder which does not break time reversal symmetry 



can not localize a single topologically protected surface 
state. Away from this energy region, there is substan- 
tial scattering caused by the step and the transmission 
drops down from the value in its absence. The total 
transmission averaged over all fc^-points is reduced for 
all energies, and one can expect that as the amount of 
disorder increases, the total transmission will be domi- 
nated by small k± contributions. We believe that these 
findings would provide a strong motivation for study of 
Sb surface with steps along T — K direction. 

In conclusion, we have performed an ab-initio study of 
the topological surface states on the Sb(lll) surface and 
their response to the presence of single bilayer high steps, 
showing excellent agreement with experimental observa- 
tions. We have identified the various scattering processes 
possible and formulated general conditions that lead to 
the formation of the dominant scattering features. This 
enabled us to confirm the fascinating chiral spin texture 
of the surface states. Resonant tunneling transmission 
across surface barriers, indicative of the extended nature 
of these states, was found. We identified phase shifts in 
the scattered PDOS at quantum well states energies and 
evaluated their life-times. The results demonstrate that 
it is possible to fully characterize the scattering proper- 
ties of the barriers with first principles calculations. Fi- 
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FIG. 3: (Color online) The allowed scattering wave-vectors 
obtained from the Fourier transform of PDOSfc^ data at (a) 
k ± b = 0, (b) k ± b = 0.16tt and (c) k ± b = 0.50tt. The panels 
on the right show the corresponding band structures at the 
same k± for the 12 bilayer slab electrodes. 
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FIG. 4: (Color online) (a) PDOS at kj_ = for surface atoms 
with a single step adjacent to a fiat region, with the step ex- 
tending along the x direction in Fig. GJb). There is an energy 
region from -60 to 20 meV with no scattering and hence no 
standing wave states, (b) Transmission at k± — 0, indicat- 
ing the perfect transmission around Ef even in presence of a 
surface barrier, the defining feature of a topological state. To- 
tal transmission also shows minimal scattering in that energy 
range. 



nally, we have shown that one can have minimal scatter- 
ing along other high-symmetry directions even in pres- 
ence of strong surface perturbations, which provides a 
unique signature for the topologically protected nature 
of these states. We believe that this can be readily tested 
in future experiments. 
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